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APPARATUS FOR REMOTE MONITORING OF A FIELD OF VIEW 



Field of the Invention 

5 The present invention relates to an apparatus for remote 
monitoring a field of view. The apparatus has particular 
application to remote monitoring of sulphur dioxide , 
volcanic ash and wind-blown dust. 

10 Background to the Invention 

There are a number of adverse atmospheric conditions that 
it would be desirable to detect. These include volcanic 
ash, toxic grases such as sulphur dioxide gas and 
15 wind-blown dust. 

Volcanic ash is a hazard to jet aircraft, causing engines 
to stall when ingested, scouring windows and the leading 
edges of the wings and causing instrument malfunctions. 

2 0 Damage to aircraft can be counted in the millions of 

dollars. Most serious aircraft encounters with ash clouds 
have been at cruise altitudes, but there is also a hazard 
to aircraft at airports affected by volcanic ash. These 
a.irports are usually close to an active volcano 
25 (e.g. Anchorage and Kagoshima) but they can also be at 
some distance from the source of the eruption due to 
atmospheric transport that brings ash into the region. 

The cost of ash hazards to airport operations is not 

3 0 known, but must be significant if the costs include those 

due to delays to landings and take-offs as well as 
re-routing costs incurred by airline operators. Currently 
there are no regulatory requirements for airport operators 
to provide warnings of ash hazards. Warnings are issued 
35 based on information from volcano observatories, 

meteorological advisories and, in some cases, radar 
observations of eruption columns. Radar information is 
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generally only reliable at the start of an eruption when 
the ash cloud is thick and usually such information is 
only available at airports in close proximity to an 
erupting volcano. For airports distant from the source of 
5 ash there are few direct observations available. Some 

observations come from satellite systems and other sources 
of information come from trajectory forecasts based on 
wind data and cloud height information. Much of this 
information is sporadic and untimely and there is a need 
10 for better detection systems. 

Other adverse atmospheric conditions include the toxic 
gases emitted by volcanoes and industrial plants. Of 
particular importance and abundance is sulphur dioxide 

15 gas. This gas is colourless, but has a characteristic 
pungent odour. Eye irritation and inflammation of th.e 
respiratory tract occurs in relatively low concentrations. 
Amounts of 6-12 ppm will cause immediate irritation of the 
nose and throat. Long term exposure can exacerbate asthma 

2 0 and can be dangerous to persons with preexisting 

cardiopulminary diseases. Thus monitoring near to strong 
sources of S0 2 (e.g. from industrial sources and at 
volcanoes) is important as is longer term monitoring at 
some distance from the source. Furthermore/ S0 2 clouds 

2 5 from volcanoes will react with water vapour in the 

atmosphere to produce sulphuric acid which can damage 
aircraft. Accordingly/ it would be desirable to be able 
to warn aircraft of S0 2 clouds. 

3 0 Wind-blown dust from desert regions or semiarid lands can 

be a hazard to aircraft/ reduces visibility significantly 
and can cause eye and throat irritation to humans. Large 
parts of the habitable earth are prone to dust storms , 
including northern Africa, the Mediterranean islands, 
3 5 southern Italy, Spain and France, southwestern USA, 

central and southern Australia, western parts of Soutli 
America, central China, Japan and south and north Korea 
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and the central deserts of Asia. The wind-blovm dust can 
also be transported long distances - dust from China has 
been detected in North America. The dust consists of 
nearly spherical particles of Si0 2 in concentrations that 
5 can limit visibility to a few 10 's of metres. 

Accordingly, wind-blown dust can be a significant hazard 
to aircraft, vehicles and the like. 

Accordingly, it would be desirable to provide an apparatus 
10 capable of providing warning of one or more adverse 
atmospheric conditions . 

Summary of Invention 

15 The invention provides apparatus for remote monitoring of 
a field of view comprising: 

a monitoring station; and 

remote sensing equipment at a site remote from 
said monitoring station, said remote sensing equipment 

2 0 compr i s ing : 

an infrared detection apparatus that monitors 
said field of view for at least two wavelengths of 
infrared radiation corresponding to an adverse atmospheric 
condition and produces temperature information based on 
25 the monitored inf raxed radiation; 

a processor for processing said temperature 
information to determine whether an alarm condition for 
said adverse atmospheric condition is met; and 

communication means for sending data to said 

3 0 monitoring station if said alarm condition is met. 

The invention provides a method of monitoring a field of 
view comprising: 

monitoring, at a remote location, a field of view 
35 for at least two wavelengths of infrared radiation 

corresponding to an adverse atmospheric condition and 
producing temperature information; 
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processing said temperature information to 
determine whether an alarm condition for said adverse 
condition is met; and 

transmitting data to a monitoring station if said 
5 alarm condition is met. 

Brief Description of the Drawings 

A preferred embodiment of the invention will now be 
10 described with reference to the accompanying drawings in 
which: 

Figure la is a block diagram of apparatus for 
remote monitoring of a field of view of a preferred 
embodiment ; 

15 Figure lb is a schematic diagram showing 

apparatus for remote monitoring of a plurality of field of 
views; 

Figure 1c is a schematic diagram which shows the 
camera portion of the apparatus of Figure la; 

2 0 Figure 2 shows the filter functions of the 

apparatus ; 

Figure 3 shows the variation of elevation angle 
with temperature difference; 

Figures 4a and 4b show variation of elevation 
25 angle with temperature difference for two different 
channel differences; 

Figure 5 shows the apparatus of the invention 
viewing a sulphur dioxide plume; 

Figure 6 shows an image produced of S0 2 plume 

3 0 using the apparatus of the present invention; 

Figure 7 shows a view of a plume of S0 2 with an 

explosion; 

Figure 8 shows an image of an explosion produced 
using an ash algorithm; 
35 Figure 9 shows a different image for the sky when 

there is no ash or sulphur dioxide present; 

Figure 10 is a temperature histogram for the 
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image of Figure 9; 

Figure 11 illustrates the Gaussian and 
thresholding technique for setting an ash alarm; 

Figures 12a and 12b show a raw and a calibrated 
5 image respectively, illustrating the dramatic effect 
calibration has on the identification of features; 

Figure 13 is a map of test sites; 

Figure 14 shows the proximity of Tavurvur volcano 
to the airport; 

10 Figure 15 has ash and visible images taken from 

Rababa; 

Figure 16 shows the alarm histogram for Figure 

15; 

Figure 17 has ash and visible images from 

15 Matupit; 

Figure 18 shows the alarm histogram for Figure 

17; 

Figure 19 is a view of an eruption from Rabaul 
Volcanological Observatory; 
2 0 Figure 20 is a graph of an alarm time- series; 

Figure 21 has ash and visible images from Hamamas 
Hotel; and 

Figure 22 shows the alarm histogram for Figure 

22. 



25 



Description of the Preferred Embodiment 



The preferred embodiment provides an apparatus for remote 
monitoring of infrared radiation from a field of view in 
3 0 order to detect an adverse atmospheric condition such as 
volcanic ash, sulphur dioxide or wind-blown dust. 



Referring to Figure la, there is shown a block diagram of 
apparatus for remote monitoring of a field of view 100. 
35 The apparatus has a remote sensing equipment 130 that has 
an infrared detection apparatus 110 and a processor 140 
provided in the form of a computer running software to 
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process data received from the infrared detection 
apparatus 110 relating to the temperature of the field of 
view. Hereafter, "temperature inf ormation" . 

5 Typically the remote sensing equipment will be located at 
a site remote from where monitoring of the atmospheric 
condition will take place. For example at a site where 
the field of view contains a volcano, whereas the 
monitoring system 150 will be at a location where it is 
10 necessary to take action in response to there being an 

adverse atmospheric condition. The communications means 
160 is typically a satellite modem that transmits the data 
to the central monitoring system. 

15 The processor 180 normally processes the temperature 

information to form a temperature difference image that 
can be used to observe the atmospheric condition and 
determine whether action needs to be taken in respect of 
the adverse conditions. 

20 

However, the images themselves are taken rapidly and 
contain a lot of information as a temperature difference 
is calculated for each pixel of the infrared detection 
apparatus 110 . Accordingly, there are potential bandwidth 

2 5 problems in transmitting the images constantly whereas it 

is only necessary for the images to be reviewed when there 
is actually a potential atmospheric hazard. Accordingly, 
the processor 14 0 processes the temperature information to 
determine whether an alarm condition is met and sends data 

3 0 to the central monitoring system 150 only when the 

condition is met. The alarm condition can be set at a 
number of levels. For example, if it is certain that 
there is an adverse atmospheric condition present or if 
interpretation of the image might lead to a conclusion 
3 5 that there is an adverse atmospheric condition present. A 
further advantage is that the temperature difference 
images do not have to be manually reviewed until an alarm 
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condition is met. 

The data sent to the monitoring station can take a number 
of forms, however, typically the data sent for further 
5 evaluation is the temperature difference image which 

caused the alarm condition to be met However, in other 
modes of operation the alarm can be a warning signal that 
allows the persons monitoring the central monitoring 
system 150 to take further action. That further action 
10 may include sending a request to the remote sensing 
equipment 130 to request a series of images. 

Figure lb illustrates a typical arrangement involving a 
plurality of remote sensing equipments 130a, 130b, 130c . 

15 Each set of remote sensing equipment 13 0 communicates with 
satellite 520 to transmit data and alarm conditions 
relating to the remote locations of which they are 
located. Satellite 520 downloads the information to an 
Internet service provider 530 that transmits the 

20 information via the Internet 532 to a second Internet 

service provider 534. The monitoring station 150 obtains 
the data from ISP 534 and can communicate via the reverse 
path with each remote sensing equipment 130. 

25 The infrared detection apparatus 110 is provided in form 
of an infrared camera 1 as shown in Figure lc. 

The camera 1 itself has a filter wheel housing 2 that has 
a window 3 that is transmissive in the infrared, a shutter 
30 4, a filter wheel cover 6 and a filter wheel. Infrared 

array housing 8 contains an infrared array 9 and a signal 
processing unit 10. 

Output from the signal processing unit 10 is on signal 
3 5 lines 11a which have an Ethernet interface to computer 140 
which can process the signals. Lines lib allow the 
shutter control signals to be received from the computer 
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140 and also temperature measurements passed to the 
computer 140. The temperature measurement corresponds to 
the temperature of the shutter 4 and is obtained by a 
contact thermometer (not shown) . 

5 

The infrared detection apparatus operates by processing 
infrared signals from the region of sky being monitored 
above at up to five pre-defined wavelengths. The 
wavelengths which are used depend on the adverse 

10 atmospheric condition which is being monitored. The 

central wavelengths and wavelength intervals are given in 
Table 1. It will be appreciated that a band of 
wavelengths surround a central wavelength/ but for 
convenience "wavelength" means the central wavelength and 

15 surrounding band unless the context implies otherwise. 
The infrared radiation measured by the camera 1 is 
linearly proportional to the resistance change in the 
detector, which is recorded and logged by the signal 
processing unit. In the preferred embodiment the infrared 

2 0 array is an uncooled microBolometer staring array of 

320 x 240 elements sensitive to radiation in the 6-14 pm 
wavelength interval is used to detect filtered radiation. 
The detection apparatus uses a filter wheel to filter 
radiation. The radiation from the sky is focussed onto 

2 5 the array by means of focussing optics in the form of lens 

5 and the field of view is a cone of up to 90 degrees. 

As indicated above, in the preferred embodiment, the 
infrared array will be an uncooled microBolometer array of 
30 dimensions at least 320 x 240 elements, but 640 x 480 

elements is also possible. There is a trade-off between 
the number of elements, cost and maximum spatial 
resolution per pixel for a fixed optical arrangement. The 
microBolometer operates on the principle that a 

3 5 temperature change produced by radiation falling on the 

detector produces a linear resistance change in the 
material. There are three types of bolometer in 



WO 2005/031321 



PCT/AU2004/001340 



commercial use: metal, semi-conductor and super-conductor. 
In the preferred embodiment either VOx and Si-based 
semi-conductor bolometers can be used as these are 
available commercially. However, it would be possible to 
5 use a cooled array with a cryogenic cooler with the 

ground-based device if the performance criteria cannot be 
met with an uncooled array. 

The camera 1 can be used to obtain temperature 
10 measurements at up to five separate wavelengths to be 
filtered and also can have a single broadband channel 
depending on what atmospheric condition is being 
monitored. The camera 1 uses a filter wheel mounted with 
circular interference filters. Table 1 provides the 
15 information for the selection of the filters. 

The filters will be 50 mm diameter germanium/ZnSe 
multi- layer interference filters mounted on a rotating 
wheel and driven by a stepper-motor. (But, smaller or 

2 0 larger diameter filters may be used depending on the field 

of view required and the focusing power of the optics) . 

The array 9 has a nominal noise temperature of no greater 
than «50 mK in the broadband channel. To achieve 
25 sufficient temperature sensitivity in the narrow band 

channels, frame averaging is employed. Table 1 shows the 
theoretical noise equivalent temperature differences 
(NEAT's) expected for various frame averaging values in 5 
narrow wavelength bands or channels and one broadband 

3 0 channel 20. Thus five narrow band channels centered 

around 7.3 Jim, 10.1 |Jm, 11 Jim and 12 pm are shown in Figure 
2 as items 21, 22, 23, 24 and 25 respectively. It will be 
appreciated by those skilled in the art that the precise 
central wavelengths and bandpasses may vary, and that the 
3 5 wavelengths and bandpasses shown in Table 1 are nominal 
and used here as indicative working wavelengths . 
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Band 


High 


Low 


NEAT 


4 


8 


16 


32 


■ 64 


7.3 


7 .05 


7.55 


1346 


673 


476 


337 


238 


168 


8.6 


8.35 


8.85 


890 


445 


315 


223 


157 


111 


10.1 


9.85 


10.35 


643 


322 


227 


161 


114 


80 


11 


10.75 


11.25 


657 


329 


232 


164 


116 


82 


12 


11.75 


12.25 


900 


450 


318 


225 


159 


113 


8-14 


8 


14 


65 


33 


23 


16 


11 


8 



Table Is Noise equivalent temperatures (NEAT, iciK) for 



different amounts of frame averaging . 

In order to produce an output indicative of the presence 
5 of an adverse atmospheric condition, it is necessary to 

calibrate or correct the raw data produced by the infrared 
array . 

The infrared camera 1 is calibrated so that the processing 
10 means 140 can produce corrected radiance values which then 
can be used to produce scene temperatures which can 
subsequently be processed using an algorithm specific to 
the atmospheric condition in order to determine the 
presence of the adverse atmospheric conditions. The 
15 calibration process consists of a pre-calibration and a 

field calibration. The calibrations correct for radiation 
from the infrared detection apparatus. The field 
calibration corrects for changes in the radiation from the 
infrared detection apparatus during operation. 

20 

Figure 12a shows a typical scene where the image is 
constructed from the raw signals, without calibration. 
Figure 12b shows the same scene after calibration and 
conversion to temperature units. Aspects of the scene not 

2 5 visible in the raw data are now clearly noticeable in the 

calibrated data. For example, roof 2 00 is now visible. 

The camera 1 provides raw digital counts as output of 
detector array 9 that have also had some corrections 

3 0 applied. These counts can be related to the scene 
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radiance through a linear calibration process, and then to 
temperature through use of the Planck function. Xt will 
be appreciated that because temperature can be derived 
from these counts, they provide temperature information* 

5 

The system uses a two point blackbody pre-calibration 
procedure that uses the output from the detector array 9 
corresponding to two cooled and heated blackbody cavities 
placed in front of the lens. The calibration equations 
10 are: 

R±,c = a±Ci, c + Jbi, (1) 
R±,h = a±C ifh + Jbi, (2) 

15 where, the subscripts c and h to the cold and hot 

blackbodies, and i refers to the channel or filter number 
being used, R = radiance, C = counts and a and b are 
coefficients corresponding to a gain and an offset 
respectively. In practice the camera views over a band of 

20 wavelengths and the response of the camera (detector and 
filter) as a function of wavelength must be known. The 
radiance is therefore related to the scene temperature, T s , 
through, 

25 R± = £ B[X, T s ]F(Z)dJl, (3) 

where X is wavelength, F is the response of the system, 
and B is the Planck function. 

30 To convert from the calibrated radiances to scene 
temperature, Equation (3) is inverted. This is a 
non- linear problem which requires a minimisation 
procedure. A series of look-up tables were generated that 
give radiances equivalent to a series of pre- specif ied 

35 temperatures. Once the measured radiance is known by 

combining (1) and (2) , the look-up table is searched and 
interpolated (if necessary) to determine the closest scene 
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temperature. The procedure is accurate to 10 mK over the 
range of observable temperatures 220 K to 330 K. 

A separate set of calibration coefficients a±, b± is 
5 developed for each pixel within the 320 x 240 image. 
Towards the edges of the image the quality of the 
calibration degrades due to image distortions and non- 
uniformity of the blackbodies. The blackbody temperature 
is measured in one place on each blackbody and non- 
10 uniformity of the temperature field will occur to some 
degree • 

A field calibration technique is used to alter the 
coefficients to account for the optics (particularly the 

15 lens) that may be heating up or cooling down and thus be 

at a different temperature to its value when calibrated in 
the laboratory. This causes an off -set in the measured 
signals. Our field calibration procedure makes use of a 
single shutter measurement just before measurements are 

20 taken at each wavelength. The shutter fills or slightly 

overfills the field of view of the instrument and provides 
a uniform radiation source to the detector. The 
temperature of the side of the shutter facing the lens is 
continuously monitored using a contact temperature probe. 

2 5 The shutter side facing the lens is blackened so that its 

infrared emissivity is high (exceeding 0.98) and uniform 
across the region 6-14 \xm. In the field, the calibration 
is performed by making a single measurement of the 
shutter, followed by a measurement of the scene and then 
30 application of the calibration equations and shutter 

measurement which accounts for the off- set generated by 
any change in temperature of the lens or other radiating 
surfaces in front of the detector. 

3 5 The calibration means also calibrates for background 

atmospheric conditions and viewing angle. That is, the 
temperature differences on a single channel will vary 
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depending on the channel measured. 

In clear and cloudy skies when there is no ash or SO2 
present, water vapour causes differential absorption of 
5 radiation in the atmospheric window between 6-14 JJin. Thus 
when comparing two channels there will be a temperature 
difference. Theoretical calculations and modelling 
studies indicate that this difference will be negative 
when the camera views the sky above the horizon. The 

10 exact value of the difference depends on the amount of 
water vapour , but also on the path length that the 
radiation traverses through the atmosphere. Figure 3 
shows the theoretical variation of the temperature 
difference (11-12 |4in) with elevation for a cloudless 

15 atmosphere containing about 3 cm of precipitable water. 
At low elevation angles the temperature difference is 
slightly negative, but gets progressively more negative 
until at around 60 degrees elevation when the difference 
decreases slowly. A consequence of this behaviour is that 

20 it is not possible to set a constant threshold for 

deciding whether infrared difference images contain ash 
affected pixels. Figure 4a shows the difference as 
determined from measurements made at Saipan. The 
variation with elevation angle mimics the theoretical 

25 behaviour. The same effect with elevation can be seen for 
8.6-12 pm temperature differences (Fig. 4b) , except that 
after 60 degrees the difference starts to increase rather 
than decrease. This is not seen in the modelling results, 
and more data are required to determine the cause of this 

30 effect. These data show more variation than the 

theoretical studies because the scene also contains clouds 
and unmodelled water vapour variations. Nevertheless, the 
temperature difference decreases with elevation angle in 
all cases studied and agrees with the theoretical 

35 behaviour. 

When the variations and temperature difference are 
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understood, temperature difference can be corrected in 
order to correct the output images so that the output 
temperature difference images correctly reflect the 
presence of the adverse atmospheric condition being 
5 detected. 

The algorithms and temperature differences which are used 
depend on the adverse atmospheric condition that is being 
monitored for in the field of view. 

10 

SQ 2 Algorithm 

Figure 5 shows a digital visible image of the camera 
viewing towards Etna volcano in the background with a 
plume of S0 2 3 0 emitted from the crater. An S0 2 image 
produced by the apparatus within 30 minutes of the digital 
image is shown in Figure 6. The colour scale 22 is drawn 
to indicate the amount of S0 2 in the plume - from yellow, 
indicating low amount, to brown, indicating high amounts. 

The background to the images comprises light areas of blue 
and green indicating a colder background whereas the 
bottom of the image which is dark in colour represents the 
ground. These areas are labelled 33 and 34 respectively. 
The left vertical axis represents elevation in degrees and 
the horizontal axis represents Azimuth. The images are 
typically produced in colour as indicated above however it 
will be appreciated that appropriate grey scale images can 
also be produced. The colour images make it easier to 
discern the plumes from the background. 

An S0 2 index is based on a 4 -channel algorithm, while the 
image itself utilizes all 5 channels, the 4 narrow band 
channels and a wide band channel. The S0 2 index is derived 
35 by: 

(1) forming the temperature difference between a 
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channel centred at 8.6 \m and a channel centred 
at lO.Ojim, label this difference as dr lr 

(2) forming the temperature difference between a 
channel at 11.0 |Jm and a channel centred at 

5 12,0 Jim, label this difference as dr 2 , 

(3) adding temperature differences <STi and dr 2 , label 
this dr 3 * 

(4) subtracting a reference value that depends on 
the viewing elevation of the camera, and has a 

10 typical range of 1-3°C to get dr 4 (this is the S0 2 

index) . 

The displayed image is produced by scaling dr 4 and 
overlaying this scaled image into the broadband image so 

15 that all pixels in the dr 4 temperature difference image 
with a scaled value in the range 1-32 are preserved and 
all pixels outside this range are replaced by the 
broadband image pixels. A suitable colour table is then 
attached to the image and a reference grid and scale are 

2 0 incorporated . 

An exemplary resulting image (Figure 6) shows S0 2 plume 30 
in yellow to brown colour, water vapour (in various 
degrees of amount) in grey colours, the background 
2 5 (colder) sky as blue and green 33 and topographic features 
34 (mountain, ground, trees etc., which are generally 
warmer than the plume) as dark grey to black. 

Volcanic Ash Algorithm 

30 

A volcanic ash algorithm which can be used with the 
apparatus may be stated as: 

dT 2 = T 1X - T 12 > AT t , 

35 

where the subscripts 11 and 12 refer to the channel 
central wavelengths and AT t is a temperature threshold that 
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depends on the water vapour content of the atmosphere and 
on the viewing elevation angle of the camera. The nominal 
value for AT t is 0°C. Data (pixels) with values above the 
threshold are regarded as volcanic ash. Data (pixels) 
5 below the threshold are regarded as not volcanic ash. 

In testing, the apparatus was also able to capture 
discrete explosions from the Stromboli Volcano. An 
example of this is shown in Figure 7. In the case of the 

10 explosion, the pyroclastic material is mostly volcanic hot 
rocks, cinders and ash and reveals itself as grey to black 
colours 41 when the S0 2 algorithm is used. An S0 2 plume 42 
is clearly shown. In contrast when the ash algorithm is 
used, that is, by taking temperature differences using two 

15 channels, specifically dr 2 , the image shown in Figure 7 is 
obtained. The colour scale now shows positive temperature 
differences in shades of orange and red, and negative 
differences as blue to yellow. In this case the algorithm 
identifies the hot rocks and cinders as positive 

2 0 differences 40 (high ash content) , and resuspended ash as 

slightly negative (similar to the material on the surface 
of the mountain slopes) . The sky has markedly negative 
differences. 

25 Wind-Blown Dust Algorithm 

Desert dust has a high silica (Si0 2 ) content and when small 
particles (diameters less than 10 pm) are suspended in the 
atmosphere they disperse infrared radiation in a similar 

3 0 fashion to volcanic ash particles. Consequently, the 

algorithm used to identify ash in the atmosphere can also 
be used to identify wind-blown dust. Dust storms are a 
frequent and global phenomenon. Most of the dust is 
confined to the boundary layer - the part of the 
3 5 atmosphere closest to the surface and generally not 

extending more than 5 km upwards. Occasionally, large 
dust storms can be transported vast horizontal distances 
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(may 1000' s of kilometres) and be lifted to heights 
greater than 5 km. Dust storms have been identified using 
passive infrared radiation from satellites. The dust 
algorithm differs from the ash algorithm in one important 
5 aspect. Since it is unlikely that wind-blown dust will 
contain any appreciable amounts of S0 2 gas (the reverse 
being true for ash) , a channel at 8.6 jam, can be used in 
conjunction with the 11 and 12 |jm channels. The dust 
algorithm thus uses three channels rather than two. The 
10 form of the algorithm is: 

9T duat = aT a . 6 + bTu + cT X2r 

where the subscripts 8.6, 11 and 12 refer to the central 
15 wavelengths (in pm) for each channel and a, Jb and c are 
constants. The nominal values of these constants are: 
a = 1, b = 1 and c = -1 . 

It will be appreciated that it will not always be 

2 0 convenient or appropriate to send data from the remote 

sensing equipment to a monitoring station. Accordingly, 
an automated algorithm can be developed in order to 
initiate an alarm. This alarm is based on analysis of the 
difference images produced in accordance with the 
25 algorithms. However, it will be appreciated that the 

images do not actually have to be produced in order for 
the alarm to be triggered, that is, the data can be 
processed instead. An example of an alarm algorithm 
suitable for volcanic ash is described, however it is to 

3 0 be appreciated that other alarm algorithms can be 

developed. 

To minimise operator intervention of the infrared image 
data and trigger the transmission of an image to the 
35 monitoring station 150, an automated algorithm has been 
developed - i.e. to transmit the image data to the 
monitoring station 150 when there is sufficient reason for 
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a user to inspect the image. The algorithm or x alarm' is 
based on a histogramming technique that takes into account 
the viewing elevation and the amount of water vapour in 
the atmosphere. In the current embodiment there are 32 0 x 
5 240 image pixels in a single difference image (ATj, AT 2 r 

AT3) . Due to cloud movement, noise, calibration errors and 
sensitivity limitations, some pixels will appear anomalous 
when there is little or no hazard within the image. The 
histogramming technique accounts for these anomalies. 

10 

In general the structure of these anomalies is very 
different to that expected from an ash cloud. However, on 
a pixel -by-pixel basis it is impossible to determine 
whether the signal is due to a camera anomaly or due to a 

15 real ash signature. Analysis of the images obtained from 
Anatahan volcano in conditions where ash was known to be 
present suggests that analyses of structure in the images 
can be used to set a threshold or alarm to indicate the 
presence of ash. To demonstrate how this can be done we 

2 0 first consider a set of images obtained in conditions 

where there was no ash or S0 2 . Figure 9 shows an image 
obtained in ash/S0 2 -free conditions viewing with an 
elevation of 20 degrees above the horizon. A colour scale 
would be present in practice on this image indicating a 

2 5 temperature range from -15 K to +10 K, with red-coloured 
pixels having the most positive temperature difference. 
To highlight the region where most ambiguity might exist, 
a grey-scale showing temperatures from -0.5 K to +0.5 K is 
included within the main colour scale. Thus grey- coloured 

30 pixels in the temperature difference image may be regarded 
as marginal, in terms of detectability . In this image 
there are some grey- coloured pixels, but the majority of 
the pixels are yellow, green to blue indicating negative 
temperature differences and hence normal conditions 

35 (i.e. clear kkies or water/ice meteorological clouds) . 
The 2 -dimensional histogram of this image is shown in 
Figure 10 . In practice the same temperature range and 
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colour scale are used for the histogram. From theoretical 
and modelling calculations we expect pixels that are ash 
contaminated to have positive differences. But, their 
actual value depends on viewing conditions, particularly 
5 the viewing elevation, and also the amount of water vapour 
in the path* Accordingly, we have determined that a 
threshold value of 0 K for ash is appropriate under most 
conditions. Figure 10 was obtained at 2 0 degrees 
elevation and as the field of view of the infrared camera 

10 is roughly 24 degrees in the vertical direction, some 

parts of the image view land surfaces. The histogram has 
prominent peaks at roughly -1 K 51 and -5 K 50 which 
correspond to clouds and clear skies, respectively. In 
this case the least negative peak has a tail that includes 

15 some positive pixels. In the corresponding image these 
pixels are due to viewing features that are low on the 
horizon and include ground targets. Such * anomalies' are 
difficult to isolate in an automated manner and could give 
rise to false alarms if a straightforward pixel 

2 0 thresholding technique were employed. 

The scheme chosen to automatically determine whether an 
image has detected ash is a statistically based method. 
This is the method of choice because by the nature of the 

2 5 problem there is often going to be a distribution of 

pixels that can be flagged as ash, within an image that 
has many pixels that are definitely ash or definitely not 
ash. In addition, because of the likelihood that pixels 
will contain mixtures, a simple threshold and binary 

3 0 decision process would be inappropriate. 

The histogram shown in Figure 10 consists of two prominent 
peaks with a spread of pixels around these peaks. If the 
detection apparatus viewed a target of constant 
35 temperature (e.g. a uniform cloud or the clear sky), then 
simply because of the fact that the camera has a wide 
field of view and there is water vapour absorption along 
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the differing paths to the target , the resulting 
difference image would be non-uniform. In practice it is 
unlikely that the sky would present a uniform target and 
even less likely that a cloud would be perfectly uniform. 
5 The combination of these effects leads to a natural spread 
in the histogram of the temperature differences, with a 
central peak corresponding to the mode temperature 
difference. For a relatively uniform scene the peak would 
be high and the spread (or standard deviation of the 
10 distribution) would be low. We have selected a Gaussian 
distribution to model the distribution. The Gaussian 
distribution in mathematical terms is, 



15 



G(Ar) = 4jexp- 



' (Ar-// Ar ) j 2 ^ 



where AT Is the temperature difference, juat is the mean 
temperature difference, Oat is the standard deviation and 
A 0 is the maximum frequency, which occurs when AT=fi A T» 
Each of the peaks (i - 2 _ a) within the frequency 

2 0 distribution (histogram plot) is assumed to be centred at 
Mat,± with a spread of cr 4r/1 . A set of Gaussian 
distributions is fitted to the frequency distribution data 
and the parameters, A 0 ,±, Mat,±, and o^i derived. The linear 
combination of these distributions is the model -fit to the 

25 data. 



The fit for the histogram data shown in Figure 10 is shown 



in Figure 11. 


Three 


Gauss ians 


were 


used 


in the fit: 


Parameter 


i=l 




i=2 




i=3 




74.2% 




24.9% 




0.9% 




-4.24 


K 


-0 .84 


K 


-0.67 K 


Gat, i 


+ 1.49 


K 


+ 0.33 


K 


+0.08 K 



3 0 The fit to the distribution although not perfect, is 
sufficient for setting the alarm for the image. A 
threshold Gaussian is set with a mean and standard 
deviation derived from modelling, and comparing this with 
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the n-Gaussian data- fit. The region between the pixels 
bounded by the threshold Gaussian mean value, and the 
overlap region between the two Gaussians (the threshold 
and the data- fit) is calculated. This area (or number of 
5 pixels) is subtracted from the number of pixels that 

exceed the threshold Gaussian mean value and lie within 
the data- fit Gaussian (see Fig. 11) . 



The ratios / 

10 / Pi-Pa 



15 



Ri = 



where P 0 ,± is the number of overlap pixels for Gaussian i, 
P± is the number of pixels that exceed the threshold mean, 
and A 0 ,± are the maxima for the Gaussian fits. The purpose 
of normalising by the maximum is to ensure that more 
2 0 weight is given to distributions that have well-defined 
and dominant peaks • 



It will be appreciated that any number of different 
statistical techniques can be used in order to determine 
25 whether sufficient numbers of the pixels relate to the 
atmospheric condition being monitored to warrant the 
generation of alarm. 



Experimental Results 

30 

During 26-30 November 2 003 CSIRO experimentally operated 
the detection apparatus at the site of an erupting volcano 
near Rabaul town, in west New Britain, PNG. The volcano, 
Tavurvur has ash- rich explosions every 10-3 0 minutes, with 
35 plumes that extend several hundred meters above the 

crater, approximately 400 m above sea level. A large 
number of infrared images were obtained at various 
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distances from the crater and employing a variety of 
viewing angles • 

The results indicate that the remote monitoring apparatus 
5 can image ash plumes and clouds and clearly discern these 
from meteorological clouds and trigger appropriate alarms. 
Results are best at closest proximity to the ash cloud, 
but good results were obtained at distances greater than 5 
km from the active crater. The ash alarm algorithm was 
10 also tested in an autonomous manner overnight from a 
distance of -8 km. 

Measurement Sites 

15 Figure 13 shows the locations of the measurement sites 
(six in all) used to image the ash- rich eruptions from 
Tavurvur. They are listed in Table 2. 



Site 


Site name 


Distance from 


label 




crater (km) 


A 


Hamamas Hotel 


3.7 


B 


Rabaul airport 
(old) 


3.0 


C 


Rababa ("hot 
springs' 7 ) 


1.8 


D 


Matupit village 


2.5 


E 


RVO 


8.0 


F 


CPL Mill 


7.4 



20 

Table 2. Site labels, names and locations used to make 
measurements of the ash- rich plumes and clouds from 
Tavurvur volcano (labeled Tavurvur in Fig. 13) . 

2 5 Tavurvur has been active since a major eruption took place 
in September 1994 which devastated the town of Rabaul and 
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destroyed the airport. A new airport, Tokua, has been 
constructed and is located about 20 km southwest of 
Tavurvur. With the crater still active, flights in and 
out of Tokua only take place in daylight hours and not at 
5 all if the winds move the ash towards the airport runway. 
Figure 14 shows a digital photograph taken from the runway 
at Tokua. A plume 60 from Tavurvur is noticeable in the 
background . 

10 The Rabaul Volcanological Observatory (RVO) operates on a 
hill overlooking the active crater and at about 8 km 
distance from it. Economic pressures in PNG have meant 
that only limited resources are available at RVO for 
operating geophysical equipment and power failures are 

15 also common . The main means of transport throughout PNG 
is by jet and light aircraft and the economy is highly 
dependent on air transportation. Thus there is an urgent 
need to monitor the volcanoes in New Britain (there are 
many) and throughout PNG. 

20 

The apparatus of the preferred embodiment operates off 
batteries for up to 16 hours and can be deployed in 
relatively hostile environments, rapidly by a single user. 
To test the ability of the instrument to distinguish ash 

25 from other meteorological clouds, the apparatus was 

deployed at a variety of locations and in a variety of 
viewing configurations. The measurements were made in an 
atmosphere with quite high water vapour amounts and at 
elevation angles varying from 10° to nearly vertical 

3 0 viewing. On many occasions the atmosphere around the 
instrument was filled with ash particles, making the 
atmosphere appear grey and causing irritation to the eyes 
and lungs . 
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Results 

( a ) Rababa 

5 The best results were obtained from Rababa , approximately 
1800 m from the crater. Note that there are also water 
clouds 70 in the image and the high frequency of activity 
has caused the atmosphere to be heavily laden with ash 
particles . 

10 

Figure 15 shows typical results obtained from Rababa (c) . 
The ash image (left-panel) correctly identifies the plume 
71 and clouds of ash from Tavurvur. Grey to black 
coloured regions of the image are identified as having no 
15 ash. The mountainside 72 is also identified as ash- this 
is not surprising since the mountain is covered in ash 
particles - 

The automatic alarm algorithm was used on all images and 

2 0 the alarm generated for the image shown in Fig. 15 is 

shown in Fig. 16. The alarm is being generated because of 
the difference between the actual histogram 210 and 
threshold histogram 211. About 43% of the pixels are 
identified as ash in the image and a clear and unambiguous 
25 ash alarm has been triggered. 

(b) Matupit village 

Good results were also obtained from Matupit village (D) 
30 about 2.5 km from the vent. The measurements from Matupit 
were made in relatively "wet" conditions and drizzle as 
well as ash fall were observed here. Results were very 
similar to those obtained from Rababa. To illustrate the 
ability of the apparatus to detect ash at high elevation 

3 5 angles, Fig. 17 shows an ash image and corresponding 

photograph of a dispersing ash cloud 75 viewed from 40° 
elevation. 
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In this case 73% of the pixels were identified as ash (see 
Fig. 18) . 

5 (c) RVO 

Ash images from RVO (E) were the most difficult to obtain 
because of the distance from the active vent (-8 km) and 
because the Observatory is perched on a hill; thus the 

10 camera could only view at relatively small elevation 

angles (-10° or less) . The combination of the greater 
distance and small elevation angles means that a 
considerable water vapour path is traversed (the water 
vapour absorption masks the positive temperature 

15 differences expected from the ash signal) . The frequency 
of eruption was so high that the air between the camera 
and the eruption column was often filled with fine ash 
particles. This has the effect of making the atmosphere 
appear grey and also causes large absorption in the 

2 0 infrared. The camera was operated continuously overnight 
at RVO. An example of the scene viewed by the apparatus 
from RVO is shown in Figure 19. Figure 20 shows a time- 
series of alarms detected by the apparatus from RVO. The 
series of triangles and circles 81 represent alternate 5- 

2 5 min sampled data and their similarity gives confidence in 

the results. 

The threshold for the alarm was arbitrarily set to a value 
of 10%. The colour threshold can be amended to suit the 

3 0 viewing conditions. The plot suggests that there were 

continuous ash emissions during the night - in agreement 
with what was observed during the day. While the highest 
alarm percentages never exceed 35%, this is a function of 
the viewing attitude of the instrument. If the instrument 
3 5 were sited closer to the volcano, then more of the ash 

would fill the field of view of the instrument and there 
would be a higher percentage for the alarm. 
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Because the atmosphere around Rabaul was constantly 
affected by ash it was difficult to obtain images which 
showed no ash. Some data were acquired looking away from 
5 the vent at meteorological clouds during the daytime which 
suggests the instrument was working as expected. Figure 
21 shows one of the images. In this image the 
meteorological cloud is mostly light yellow or blue-green, 
suggesting no ash. The pixels coloured yellow in the 
10 bottom right of the image correspond to very low elevation 
angles and ground targets that often give temperature 
differences slightly greater than zero. The corresponding 
alarm histogram was not triggered by these data (see Fig. 
22) . 

15 It will be apparent to a person skilled in the art that 
these and many other variations fall within the scope of 
the present invention. 



